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Abstract 

Isospin structure of gluon mediated or strong penguin is significantly altered 

when the full electroweak corrections are included. This has the consequence 

that some previous analyses which relied on a simple isospin structure in 

charmless B decays become inapplicable. We present the general Hamiltonian 

in next-to- leading order QCD, and illustrate our conclusion quantitatively for 

B — ► 7T7T and B — > Kir decays in the factorization approximation. Some 

remarks on CP asymmetries in B decays are also made. 
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Flavor changing one loop processes mediated by a gluon play an important role in un- 
derstanding the Standard Model (SM) The CP violating parameter e'/e in the K 
system is estimated by evaluating the flavor changing process s —>■ dqq, where q = u or d. 
It is by now well known that substantial corrections arise in this estimate from inclusion of 
the full electroweak effects jS],^| (the 7, Z penguins and the "box" contributions). These 
additional electroweak corrections add AI = 3/2 contribuitons to the AI = 1/2 pure gluon 
exchange, and substantially reduce the estimate for e'/e. A complete next-to-leading order 
(NLO) calculation now exists for the full electroweak contributions for such flavor changing 
process ||. In the K system the enhancement of AI = 1/2 amplitude over AI = 3/2 by a 
factor of about 20 in K — > tttt decays further amplifies the electroweak contribution to e'/e. 
One might think that in the B system electroweak corrections to gluon exchange may be 
negligible. This in fact is incorrect. The large mass of the top increases the relative contri- 
butions of the electroweak corrections and has the effect of causing significant reduction in 
the estimates for gluon mediated decays. Our recent estimates for b —> s<j) process suggests 
a reduction in the rate of 20%~ 30% ||£§]. The gluon mediated process b — > sqq where 
q = u , d , s is pure A J = 0. Inclusion of the full electroweak corrections will introduce a 
significant AI = 1 admixture. Any analysis that relies on penguin contributions being pure 
AI = is therefore suspect. We will comment on an interesting analysis assuming the above 
property which was recently published in Physical Review Letters 0. Similarly, one might 
worry about the wrong isospin admixture in b — > dqq process where the gluon exchange is 
pure AI = 1/2. Extraction of penguin effects from the tree effects are important in B — > tttt 
process J8],|9| because of the need to remove penguin contamination to make a measurement 
of the angle a in the unitarity triangle of the KM matrix. 

We shall first present the full Hamiltonian describing flavor changing processes in the 
NLO approximation. The estimate of matrix elements to quantify our results can only 
be done in the context of a model in our present calculational limitation. We shall use 
factorization approximation and some models for form factors when necessary. 

The effective Hamiltonian responsible for charmless B decays can be parametrized as 
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5ab=i = -^[VubV^Ol + c 2 O u 2 ) - V tb V t * q (HOi] + H.C. 

V 2 i=3 



(1) 



where the Wilson coefficients (WCs) q are defined at the scale of fi ~ m b , and Oj are defined 

as 



q' 
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Here q can be cf or s, and O2, 0\ are the tree level and QCD corrected operators. O^-q 
are the strong gluon induced operators, and operators O 7 _i are due to 7 and Z exchange, 
and "box" diagrams at loop level. We refer to the latter as electroweak penguin operators. 
The summation on q' in all the operators are over u , d , s quarks. The J2 q ' terms in 
Q3-6 transform as isosinglet, while the presence of e g > in Q7-10 makes the J2q' e g' terms 
transform as a mixture of / = and 1 = 1 components. Consequently, have two 
isospin components, the strong penguin has a single isospin component, and the electroweak 
penguin again contains two isospin components. 

Naively, the electroweak penguin contribuiton compared with the strong penguin con- 
tribution is suppressed by a factor a em /a s and therefore one might be tempted to neglect 
them. For large top quark mass, this is, however, no longer true because there is a term in 
the electroweak penguin contribution in which the WC is proportional to the squre of the 
top quark mass. 

The WCs Q at the scale fi = m b are obtained by first calculating the WCs at \i = mw 
and then using the renormalization group equation to evolve them down to m b - We carry 
out this analysis using the next-to- leading order QCD corrected WCs following Ref |J . We 
use in our analysis, a s {m z ) = 0.118, a em {m z ) = 1/128, m t = 174 GeV and m b = 5 GeV, 



and obtain || 

ci = -0.3125 , c 2 = 1.1502 , c 3 = 0.0174 , c 4 = -0.0373 , 
c 5 = 0.0104 , c 6 = -0.0459 , c 7 = -1.050 x 10" 5 , 

c 8 = 3.839 x 10~ 4 , c 9 = -0.0101 , c 10 = 1.959 x 10" 3 . (3) 

We see that the coefficient eg arising from electroweak penguin contribution is not much 
smaller than coefficients of the strong penguin. The major contribution to eg arises from Z 
penguin. The dependence of q on m t is given in Ref. 0. 

Our coefficients are given in a regularization independent scheme which requires that the 
matrix elements be renormalized to one loop level for consistency. These one-loop matrix 
elements can be rewritten in terms of the tree-level matrix elements < Oj >* ree of the 
effective operators, and one obtains |p|,^J10| 

< cA >= E CiOi) [<% + ^ml + ^m%] < O j >^= Y,<<Q l > tree ■ (4) 

Here m s,e are 10 x 10 matrices which we have evaluated. Expressing the effective coefficients 
d { which multiply the matrix elements < Oi >* ree in terms of Cj, we have 

c[ = Ci c' 2 = c 2 

c 3 = c 3 — P s /3 , c' 4 = C4 + P s , C5 = C5 — P s /3 , Cg = Cq + P s , 

C ? = C7 + P e , Cg = C 8 , Cg = Cg + P e , C 1Q = C10 • (5) 

The leading contributions to P s ^ e are given by: P s = (« s /87r)c2(10/9 + G(m c , /i, q 2 )) and 
P e = (a em /97r)(3ci + C2)(10/9 + G(m c , fi,q 2 )). Here m c is the charm quark mass which we 
take to be 1.35 GeV. The function G(m,fi,q 2 ) is give by 

2x t f 1 , m 2 — x(l — x)q 2 . . 

G(m,/i,g 2 ) = 4 / rr(l - a;)dxln ^ — , (6) 

and q 2 is the momentum of the exchanged particle. In factorization approximation, the final 
state phase only arises from the imaginary part of G. 
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a) B — > nit and extraction of sin2<x 

In this case q = d, and the effective Hamiltonian contains AI = 1/2 and AI = 3/2 
amplitudes. The AI =1/2 component will contribute to decay amplitudes with 1 = and 
/ = 1 in the final states, while the AI = 3/2 component will contribute to the amplitudes 
with 1 = 1 and 2 in the final states. However, Bose symmetry requires tttt to be in I = 
or / = 2 state. The KM matrix elements in front of 0± t 2, the strong penguin, and the 
electroweak penguin contributions are comparable, while the WCs of the strong penguin 
and the electroweak penguin are all smaller than the ones for 0\^- 

The extraction of a, one of three angles of the unitarity triangle defined by the KM 
matrix elements, involves the study of CP asymmetry in B — > tttt and B — > tttt channels. 
If the pegnuin contributions are neglected, this extraction is straightforward. However, if 
penguin diagrams make a significant contribution, then the interpretations of the results 
become complicated. An isospin analysis of B — > tttt has been presented by Gronau and 
London || , where a method of removing penguin contribution was provided. An important 
assumption in this method is that the AI = 3/2 amplitude involved in B — > 7r ± 7r° decays 
arises purely from the operators 0\ i. Since we now have the electroweak penguin operators 
with AI = 3/2 part, one can ask if this method is still valid. We calculate the ratio of the 
penguin to the tree amplitude in the factorization approximation, 

A tree (B + - tt+tt ) = -^v: b V ud (c[ + c' 2 )(l + OT , 



G 3 

A penguin (B + -> tt+tt ) = -J=V t * b V td -[c' 7 + £4 + 4 + £c' 10 + £4 + c ' w 



+ l(& + c' 8 )(2W + X)]T, (7) 

where £ = 1/N with N being the number of colors, W = m^/lma + m u )(mb — m^), X = 
ml/2m d (m b - m d ), and T = i/ 7r o(/+ 7r _ (m| - ml) + f Bn -m 2 n ). Here fa = 93 MeV is the 
pion form factor. For our numerical calculations we will use m u = 5.7 MeV, md = 8.7 MeV, 



and the form factors / calculated in Ref. IIJI2]]. We will treat £ as a parameter and use 



the experimental data favored number 1/2 [12 . 
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In factorization approximation, the ratio of the tree and penguin amplitudes do not 
depend on the form factors. We find 

'twB^C?' " W-I/IV-I). (8) 
Here we have neglected a small contribution about 0.5% due to strong penguin from isospin 
breaking effect for rrid ^ m u . 

We see that the assumption made in Ref. f| is quite good. We have also calculated 
the ratio of penguin to tree amplitudes for B° — > 7r + 7r~ and B° — > 7r°7r°, and find them 
to be about 7%(| Vtd|/|V^|) and 23%(|Vtd|/|V^b|), respectively. The larger penguin effect in 
B° — > 7T 7r° can be understood from the color suppression of the tree contributions to this 
process in factorization approximation. 

The method suggested in Ref. || requires measurements of all the decay amplitudes, 
B + — > 7T + 7r°, B° — > 7T + 7r - , and 5° — > 7r°7r°. Experimentally it will be very difficult to 
measure the decay amplitude for B° — > 7r°7r° to the desired accuracy for a determination 
of a. The CP asymmetry in S° — > 7r + 7r~ decay will be measured at the B factory. We 
estimate below the error in a determination from this asymmetry in presence of penguin 
contamination. The experimental measureable quatity is Im£ + _ where |TB| 

^ A(B° - 7T + 7T-) ' 1 } 

and e~ 2l/3 = V t * b Vtd/V t bV t * d . The 5° decay amplitude in general has the form, e l7 a tree + 
e~ tf3 e t5 a penguin . Here at ree and a penguin are the tree and the penguin contributions, respec- 
tively, 7 = Arg(V* b V us ), and 5 is the difference between the tree and penguin rescattering 
strong phases. In the absence of penguin contributions, Im(£ + _) = — sin2o;. Including the 
penguin contributions, we find the deviation from — sin2a to be 

A(sin2o;) = Im(£+_) + sin2a 

2cos<5sin3a — 2cos(<5 + a)sin2a — i?sin2a — i?sin4a , . 

1 + R 2 - 2ifcos(5 + a) ' ( } 

where R = a penguin / a tree . For B° — > 7r + 7r~ decay, R is about 7%. We find that 

|A(sin2a)| < 0.14 for all values of a and 5. The error in a determination could be quite 



large. For example for a = 45° this error leads to a determination of a off by 13°. 



b) Isospin analysis in B — ► Kn decays. 

We now consider the general Hamiltonian with q = s. The operators Oi t2 transform as 
AI = and 1, the strong penguin operators 3 _q transform as AI = and the electroweak 
penguin transform as AI = and 1. The operators 0\ and 2 are suppressed by a factors 
iKTfc^sl/l V^Vts| ~ 1/50, and penguin operators dominate in spite of smaller WCs. It was 
show recently that isospin analysis of B — > Kit can determine the angle 7 in the KM 
unitarity triangle |7j]. The idea proposed in Ref. f?j] was to isolate AI = 1 transition in 
B — > Kir decays, and if one could assume that this amplitude arose from the operator 0\ 2 i 
one would know its weak phase. The presence of electroweak penguins, whose contribution 
to this process is comparable to 0\ 2) makes this idea entirely unworkable. We shall set up 
the general analysis of the problem with the inclusion of electroweak penguin, and arrive at 
quantitative estimates based on factorization approximation. 

The decay amplitudes for the charged B decay into a kaon and a pion have the following 
isospin decomposition 

A{B+ - K\+) = ^A 3/2 - ^A 1/2 , 

A{B+ - K+n°) = ^A 3/2 + ^A 1/2 , (11) 

where A4 indicates the amplitude with / = i in the final state. 
Consider the combination with pure 7 = 3/2 final state, 

A(B + -> K\ + ) + V2A(B+ -> K + n°) = V3A 3/2 . (12) 

In general we expect Ai to have the structure in Wolfenstein parametrization, 

A-i ~^y^\^ub^'iis\\^'tree,i^ & ' ^penguin, i& ' ] 5 (13) 

where ^T(P),i is the tree (penguin) strong rescattering phase. We resort to the factorization 
approximation to estimate the magnitudes of the tree and penguin contributions. We find 



V3A 1/2 = ^{V: b V us [-(c[ + &)C - (£4 + c' 2 )(T' + 3A>) 



where 



- VS,Vu[-S(^ + c' A ){T' + A') - 6(£4 + C ' 6 )(YT' - ZA>) 

+ ~(4 + &)C + 6(£c' 7 + c' 8 )ZA' - 3(^4 + c' l0 )A' ~ |(4 + £ c io)^']} , 

f^3/ 2 = ^{ou-(4 + ecyc - + 4)t'] - v: b v ts [\{c' 7 + e4)c" 

- 3(e4 + 4)ft' - + - f (4 + £4o)c"]} , (i4) 



(m b - m u ){m s + m u ) ' (m b + m u ){m s - m u ) 
T' = ifK{ftA m B ~ m l) + fB^m 2 K ) , 



4), 



C = ifAfBR^B ~ m )<) + JbK™\ 

A' = -if B (ftA™K - <) - f^ml) . (15) 
Note that if penguin contributions are neglected, we would find in SU(3) limit, 



\M,2 = ^-A(B+ - vr+vr ) , (16) 



as assumed by Ref . |?J • However, we now find 

A 1/2 = ^f\V* b V us \(-0.7A7e^e fST ^ 2 + 7.305e^ 1 / 2 )(GeV 3 ) ; 
V 2 

= -^|y u *y us |(-l.055e i7 e^ 3 / 2 + 0.843e^ 3 / 2 )(GeV 3 ) . (17) 
v2 

In the above we have used: V us = 0.221, \V ts \ = \V cb \, \V ub /V cb \ = 0.08, f K = 158 MeV, and 
m s = 0.175 GeV, fg = 200 MeV, and the form factors calculated in Ref. JTT|. The first term 



is the contribution from the operators Oi^, and the second term is due to the electroweak 



penguin. We have also carried out a calculation using the from factors in Ref. |12[]. We 
find that the ratios of the tree and penguin amplitudes are about the same as the previous 
ones although the magnitudes are about 30% larger in the latter case. It is clear that the 
electroweak penguin contribution to A 3 / 2 is comparable to the contribution from 0\ i- This 
makes the analysis of Ref. [0] invalid. 
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We would like to remark that B + — > K + tt° may still be a good place to look for CP 
violation [|Kj]. The CP asymmetry defined by the following equation, 

_ r(B+ -> Kn) - T(B- -> /fir) 

can be estimated from eq.(14). The quark level strong rescattering phases is zero, 5p,i/2 
and 5p )3 / 2 are estimated to be 13° and —0.23° respectively, using q 2 = ml/2 in eq.(6) with 
either from factors in Ref.ll or Ref. ||12|| . We find the asymmetry A asy can be as large as 
10% depending on the angle 7. This estimate relies on the quark level strong rescattering 
phases. There may be large long distance contributions to the strong rescattering phases 
from the final state hadron interaction. These contributions may dilute the CP asymmetry 



estimated here |14|,|15| . Unfortunately one does not know how to calculate these long distance 
contributions. 
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